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Abstract 

We study electron transport through a triple quantum dot in ring or interferometer configuration. In particular, 
we analyze the influence of a gate voltage that detunes one of the dots, such that it becomes off- resonant. In this 
i regime, interference effects fade away, i.e., the current becomes independent of a penetrating flux. Despite the absence 
^ of interference effects, the off-resonant dot causes intriguing noise properties which we characterize by the full-counting 
*™~ 5 statistics of the transported electrons. It turns out that the detuning causes strong electron bunching. Analytical results 
for limiting cases support this picture. A possible application is the construction of current sources with widely tunable 
,— ' noise properties. 

1_ i ' Keywords: quantum dots, quantum transport, full counting statistics 
^3 PACS: 73.23.-b, 05.60.Gg, 74.50. +r 



CZ3 



C3 



I 

C 

O 
o 



(N 
> 
(N 

in 

m 

(N 
O 
O 



X 



1. Introduction 

During the last decade, a huge effort has been made to 
understand the conduction properties of quantum systems 
that consist of only a few discrete levels. Spurred by Feyn- 
man's vision of "plenty of room at the bottom" [![, it has 
e.g. been proposed to use single molecules as elements of 
future electronic circuits [2J . Although this task is far from 
being accomplished, molecular electronics already became 
jan established field d 0| • Apart from their technological 
promises, conducting molecules may also serve as tools for 
the implementation of fundamental physical phenomena. 
One example is the intriguing effect of tunnel suppression 
by the purely coherent influence of an ac field [5j, which 
is stable even in the presence of Coulomb repulsion @. 
Coherent destruction of tunneling leaves its fingerprints 
in the transport characteristics of laser-driven molecular 
wires, where an ac field may suppress the current and its 
fluctuations 0, S, H E3] • 

From a theoretical perspective, molecular wires share 
many features with quantum dots. In particular, the elec- 
tronic structure of both systems consists of discrete states. 
Owing to this fact, coherently coupled quantum dots may 
be considered as "artificial molecules", despite the fact 
that their energy scales are several orders of magnitude 
smaller than those of real molecules. For roughly one 
decade, the state of the art has been to couple just two 
quantum dots coherently [Tl|, [H, EH, while trip le quan- 



tum dots have been realized only recently [14j, [15|, |16 |. 
Both double quantum dots [17| and triple quantum dots 
can be constructed such that electrons coming from 
the source may proceed on two different paths towards the 
drain. There they interfere constructively or destructively, 
depending on the setup and a possible flux enclosed by 



the interfering paths. Moreover, there exist dark states 
which are superpositions of states that are decoupled from 
the drain and, thus, block the electron transport flil ll^]. 
Such blockade may be resolved by excitations with proper 
ac fields [H [2l[ which also create spin correlations be- 
tween transported electrons [22 1. 

Here we investigate the conduction properties of a triple 
quantum dot in ring configuration as sketched in Fig.[T] A 
penetrating flux allows the operation as Aharonov-Bohm 
interferometer. We focus on the influence of a gate voltage 
that allows shifting of the levels of dot 2 out of resonance. 
The intuitive expectation is that for strong detuning, the 
two other dots govern the transport process, such that the 
triple quantum dot behaves like a double quantum dot. 
While this is indeed the case for the current, electrons 
temporarily trapped in dot 2 may cause Coulomb block- 
ade and, thus, interrupt the transport such that the elec- 
tron flow becomes avalanche-like. This is reflected in the 
zero-frequency component of the current-current correla- 
tion function which significantly exceeds the value for a 
Poisson process. With this criterion, avalanche-like trans- 
port has been predicted for quantum dots coupled to a 
harmonic mode [23l . 24 [, multi- level quantum dots (2S| . 
triple quantum dots in the Kondo regime [261 ] . and also 
when destructive interference effects strongly suppress the 
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current 

An established theoretical tool for characterizing such 
correlated transport is full counting statistics H 31, 32| 
which provides the complete information about the dis- 
tribution of the transported charge in terms of the corre- 
sponding cumulants. For a Markovian master equation as 
employed below, it is possible to express these cumulants 
as derivatives of a particular eigenvalue of the Liouville 
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Figure 1: Triple quantum dot in ring configuration coupled to elec- 
tron source and drain. The onsite energies of dots 1 and 3 are 
ei 3 = 0, while the onsite energy of dot 2 can be tuned by a gate 
voltage such that £2 = eV ga tc- The bias voltage V is assumed to shift 
the chemical potentials symmetrically, such that (i^ = eV/2 = — jUr 

operator augmented by a counting variable [32]. Gener- 
ally, this requires the computation of derivatives of high 
order. As soon as one has to rely on a numerical treat- 
ment, practical calculations may represent a formidable 
task. Recently, Flindt et al. have found a way to circum- 
vent this difficulty [33|,[34j]. Based on Rayleigh-Schrodinger 
perturbation theory, they derived a scheme that allows one 
to recursively compute the full counting statistics. We use 
this approach for our numerical solution. 

Our paper is organized as follows. In Section [5] we in- 
troduce a system-lead Hamiltonian and derive a master 
equation formalism for the computation of the full count- 
ing statistics. The numerical results for the current and the 
zero-frequency noise, presented in Scction[3J provide infor- 
mation about the super-Poissonian nature of the transport 
process. In Section 2J we derive for the cumulants up to 
fourth order analytical expressions valid in the low-bias 
regime. 

2. Model and master equation approach 

2.1. Dot-lead Hamiltonian 

We consider a triple quantum dot with a ring-shaped 
geometry as depicted in Fig. [TJ which is described by the 
Hamiltonian 

#TQD =^2eic\ci + T Y1 ( C * C J + °3 C ^ + U 12 C \ C A C ^ 

(i) 

where i,j = 1,2,3 label the dots, while c\ and Cj denote 
the corresponding creation and annihilation operator of a 
spinless electron. The first term refers to the onsite en- 
ergy €i of an electron on dot i, for which we assume that 
ei = £3 = 0, while ti = eV &Sj tc can be tuned by a gate 
voltage. The second term describes electron tunneling be- 
tween the dots, for which we assume that the tunnel ma- 
trix element T is the same for all three possible transitions. 
The last contribution models inter-dot Coulomb repulsion 
with strength U . In the present case, this term can be 
written in the form ^UN(N — 1), where N denotes the 



total number of electrons on the three dots. We assume 
that the intra-dot interaction is practically infinite, which 
corresponds to modeling each dot by a single level. 

In a ring-shaped geometry, the electrons possess two 
paths for traveling from the source to the drain. Thus, in- 
terference plays a role. This interference can be controlled 
by a magnetic flux $ penetrating the ring. We capture 
this by the replacement c\cs —> c\c% exp(27r$/$ ), where 
$0 = h/e denotes the flux quantum. Note that since we 
consider only spinless electrons, the effect of an associ- 
ated Zeeman splitting is ignored. Nevertheless, the flux 
dependence of the current provides information about the 
relevance of interference effects. 

Dots 1 and 3 are attached to electron reservoirs with 
chemical potentials /iL and /iR, such that the bias voltage 
fulfills eV = /iL — jttR. They are described by the lead 
Hamiltonian 

-ffleads = ^ e q c \q c £qi (2) 
t,q 

where I = L,R, and the expectation value {c\ q ci' q >) = 
f(e q - Ht)5a>5qq> = fe(e q )SwS qq > with chemical potential 
Ht and the Fermi function f(x) — [exp(x/fcsT) + l] _1 . The 
dot-lead contact is established by the tunnel Hamiltonian 

tftun = ^(^Lgc[ gCl + VkgC^Cs) + h.C. (3) 
1 

We assume within a wide-band limit that all effective cou- 
pling strengths I^(e) = 2TrJ2 q |V« g | 2 <5(e - e q ) are energy 
independent and that the setup is symmetric such that 

r L = r R = r. 

2.2. Markovian master equation 

The derivation of a master equation starts from the 
Liouville-von Neumann equation for the full density oper- 
ator, ihR = [i?TQD + -Stun + -ffieads, R]- By standard tech- 
niques [35j |. one obtains for the reduced density operator 
of the central system, p = tri ca( j s i?, within second-order 
perturbation theory the Bloch-Redfield equation 

i 1 f 00 

P = --[H T qd,p] - triads / dr[i? t un, [H tU n(-r) , R}} 

= Cp, (4) 

which can be evaluated under the factorization assumption 
R = Pleads, ® P- The tilde denotes the interaction pic- 
ture with respect to the uncoupled Hamiltonian, X(t) = 
U 1 (t)XU (t), where U (t) = exp{-i(#TQD + H leads )t/h}. 
We proceed by inserting the dot-lead tunnel Hamiltonian 
([3]) and evaluate the trace over the lead states. In order 
to cope with the interaction picture, we decompose the re- 
sulting master equation into the eigenstates \a) of iJxQD, 
where E a denotes the corresponding eigenenergy. This al- 
lows us to evaluate the r-integration, which yields a delta 
function and a principal value term. Neglecting the latter, 



we obtain for the density matrix elements the equation of 
motion 



Pa/3 



(5) 



a',/3' 



with the incoherent dot-lead tunneling given by 
= E ^{{a\c\W){P'\ci\P)h{E a -E al ) 

4=1,3 

+{a\c\\a'){p'\c i \(3)f t {Ep-E f) ,) 
-Y,{P'\ c Al){l\4\P)h{E 1 -E^)5 aal (6) 

7 

-^(a|c,| 7 )(7l^|a')^(^7-^')^} 

7 

+ (c£,4,fe) (c],q, 1 - /«), 

where the leads now are labeled by the number I of the 
dot to which they are attached. Electron tunneling from 
the leads to the dots is described by the explicitly written 
terms, while the replacement in the last line yields the 
terms for tunneling from the dots to the leads. 

The solution of the master equation ([5]) contains the 
full information about the state of the central system and 
provides all corresponding expectation values. However, 
we are interested in the statistics of the charge transported 
after the initial preparation, which is an expectation value 
of lead operators. Therefore we have to generalize the mas- 
ter equation formalism by introducing a counting variable 
which allows one to keep track of this information. 

2.3. Full- counting statistics 

The counting variable x is defined via the moment gen- 
erating function 



4>{ X ,t) = (eaq>(ixiV a ))t, 



(7) 



where = c l q c 3q 1S the electron number operator 
of the right lead, while the angular brackets refer to the 
expectation value at time t. The fcth derivative of 4>{x, t) 
with respect to i% at % = obviously is the moment (iVp) 
of the electron distribution in the right lead. The cumu- 
lants of the distribution are defined as the corresponding 
derivatives of ln</>(x,i). For a Markov process, they even- 
tually become linear in time. Their time-derivatives at 
large times, 



d d k 
dt d(ix) k 



ln#x,t) 



X=0,t— f oo 



(8) 



are the stationary current cumulants and characterize the 
transport. The first and the second cumulant, C\ and C2, 
are essentially the current / = eC\ and its zero-frequency 
noise S = e 2 C2, respectively. Their ratio, the Fano factor 
F = C2/C1, represents a dimensionless measure for the 



noise of the transport process [36| . It is defined such that 
for a Poissonian process F = 1. 

Our goal is now to find a reduced master equation 
that allows one to compute <p(x,t). We start again from 
the full density operator R, but now multiply it with the 
operator exp(ixA'ft) before tracing out the leads. This 
yields the generalized reduced density operator P(x, t) = 
trieads{exp(ixiVR)i?} which contains information about the 
electron distribution of the drain and fulfills the trace 
condition tr P — cj>. Using the commutation relations 
[N K ,A] = A and [JV R ,At] = -At, where A = £ g V q clc 3q , 
we obtain the master equation 

P( X ,t) = C x P( X ,t) (9) 

with the augmented Liouvillian 

£ x = £+{e i x-l)j+ + (e- ix -l)J- (10) 

and the particle current superoperators . After some 
algebra, we find for the latter the expression [37( 



eTp 



J~P=^ J q dr J dee-^cl(-r)pc 3 f e (e)+h.c., (11) 

and Sf + is formally obtained from J~ by the replacement 

(c3' c 3>/r) (c3,4' 1_ /r)- Tne superoperator J+ de- 
scribes tunneling from dot 3 to the right lead, while J~ 
corresponds to the reversed process. We proceed as for 
the derivation of the master equation (J5J) and decompose 
the current superoperators into the eigenstates of i?TQD 
which yields 



ftp,*? ~{a\c 3 \a'){a'\4m 



(12) 



{2 — fn(E a —E a ') — fn,{Ep—Epi)}, 



(13) 



{f R (E a -E a ,) + f R (Ep-Ep,)}. 



In the long-time limit, the dynamics of P(x, t) is gov- 
erned by the eigenvalue of C x with the largest real part, 
denoted as A(x). We assume that A(x) can be uniquely 
identified by its limit \iva. x -^o A(x) = 0, i.e., it corresponds 
to the stationary solution of the Liouvillian C — C x ^q. 
Then, P — A(x) exp[A(x)i], with A(x) be the correspond- 
ing "eigenoperator" of C x . It is straightforward to see 
that ln0(x,i) = hitr A(x) + X(x)t- In the long-time limit, 
the contribution of the normalization factor is not rele- 
vant and, thus, we can conclude that A(x) is the current 
cumulant generating function [32]. It can be written as 
the series 



SI 

A(x) = £Tr(ix)< 



fc=l 



(14) 



The remaining task is now to compute the proper eigen- 
value of C x and its derivatives with respect to x- I R niany 
cases, the reduced density matrix has a sufficiently small 



3 



dimension or possesses symmetries, such that one can con- 
tinue with analytical calculations. As soon as one has to 
resort to a numerical treatment, however, one faces the dif- 
ficulty of numerically computing derivatives. This can be 
avoided with the recursive scheme developed in Ref. [33| 
even for non-Markovian master equations. Here we restrict 
ourselves to the Markovian limit. 

Since we are interested in the derivatives of the cumu- 
lant generating function at x = 0, we can treat x as small 
parameter and employ perturbation theory. Then the se- 
ries (fT4"]l for the eigenvalue A(x) corresponds to the usual 
ansatz for the eigenvalue. It can be computed by the re- 
cursion derived in Appendix |Appendix A see Eqs. (|A.6[) 
and (|A.7p . Upon setting in these expressions Ek — Ck/k\, 
V k = [J+ + (-l) k J-]/kl, and P k = \</> k )/k\, we find 



k-l 



k'=0 
k-l 



(15) 



k'=0 



(16) 



The iteration starts from the stationary solution of the Li- 
ouvillian, Pq = |0o) = Poo, which implies Eq = Cq = 0. 
Multiplication with the corresponding left eigenvector (cj>o \ 
must correspond to computing the trace. This becomes 
clear when one notes that the Liouvillian is trace conserv- 
ing, i.e. ti £X — for any operator X. Consequently, 
Q = (1 — poo tr) is the projector on the subspace perpen- 
dicular to poo. In this subspace, the Liouvillian possesses 
the pseudo-inverse Q/C. Applying it to any operator, i.e. 
computing (Q/C)B = X, is equivalent to solving the lin- 
ear equation CX = B under the constraint tr A = 0. 

In particular, we obtain for the first two cumulants, 
which are the current and the zero-frequency noise, the 
known expressions 



23] 



(17) 



(18) 



I = etr{J+- J-) Poo 
S = e 2 ti(J+ +J-) Poo 

- 2e 2 tr(J+ -J-)j{J + - J- )p c 
3. Numerical observations 



An intuitive picture for the transport properties as a 
function of the gate voltage can be provided for the lim- 
iting cases in which the detuning of dot 2 is either much 
larger or much smaller than the inter-dot tunneling: In 
the limit leV^td <C T, all three dots are in resonance and, 
thus, the electrons may take with similar probability two 
different routes, namely |1) — > |3) and |1) — > |2) — > |3). 
Therefore one expects interference to be relevant, such that 
the current can be modified by a magnetic flux penetrat- 
ing the ring. If the gate voltage is large, by contrast, i.e. 
for |elgato| 3> T, dot 2 is off-resonant and, thus, should be 
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Figure 2: Current as a function of the gate voltage V ga t c = £2/e for 
different magnetic fluxes <E> and bias voltages. The inter-dot tunnel 
coupling and interaction are T = 2T and U = 80r, respectively. 
Inset: Corresponding visibility for eV = 120f\ 



of minor relevance. Therefore, the transport properties of 
the setup are expected to be essentially those of a double 
quantum dot formed by dots 1 and 3. Our numerical re- 
sults, however, will demonstrate that this is only true for 
the average current, while the current noise is significantly 
altered by the presence of dot 2 even when far detuned. 

As a first step, we solve the master equation numeri- 
cally. In doing so, we identify parameter regimes in which 
the stationary state of the triple dot is dominated by a few 
eigenstates. This allows us to reduce the complexity of the 
master equation, such that we can achieve an analytical 
treatment that provides insight to the transport mecha- 
nism. For the numerical solution, we use the dot-lead cou- 
pling r as energy unit. For the typical value T = 10/ieV, 
the corresponding current unit reads eT/h = 2.43nA. 

3.1. Stationary current and occupation probabilities 

The stationary current as a function of the gate voltage 
^gato = £2/e is shown in Fig. [3J It demonstrates that for 
Nl ^ T, the current indeed is significantly flux dependent 
as conjectured above. For flux $ = 0, the current assumes 
its minimum. This can be explained by the fact that for 
$ = = 62, the state (i/'dark) = (|1) — |2))/v2 is an eigen- 
state of the triple dot Hamiltonian and is orthogonal to 
state 1 3). Therefore it is not directly coupled to the right 
lead, which has the consequence that IV'dark) may trap an 
electron which henceforth blocks transport 38]. This is 
formally related to an atomic state that is decoupled from 
a light field and, thus, remains invisible (39j |. 

As a criterion for the relevance of interference, we com- 
pute the visibility 



where /„ 



(19) 



max$ is the maximal current upon 



flux variation for a given gate voltage, and I m i n is defined 
accordingly. The inset of Fig. [5] demonstrates that for 
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Figure 3: (a) Current as a function of the bias voltage eV = fi^ — /ir 
for various interaction strengths. The inter-dot tunnel coupling and 
the gate voltage are T = 2T and V ga tc = £2/e = 40r/e, respectively, 
(b) Probability that for U = 80r, the asymptotic state of the triple 
quantum dot contains N electrons. 



1^2 1 3> T, interference plays a minor role. The reason 
for this is that when dot 2 is strongly detuned, transport 
through this off-resonant dot requires co-tunneling and, 
thus, the path via dot 2 has a significantly lower probabil- 
ity than the direct path. 

Henceforth we focus on the region without interference 
effects. There a main feature of the current is that it ex- 
hibits plateaus, see Fig. EJa). This is consistent with the 
usual Coulomb blockade scenario in which the bias and the 
gate voltage, determine the maximal number of electrons 
that can reside on the dots. Accordingly, the probabil- 
ity for having a particular dot occupation changes at the 
steps, as can be appreciated from Fig. &h). Moreover, 
an electron can tunnel from a left lead to the dots only if 
it has sufficient energy to compensate the Coulomb repul- 
sion of the electrons that are already in the triple quantum 
dot. Thus the occupation with N electrons (N = 1,2,3) 
is possible only if one chemical potential is larger than 
(N — 1)Z7. For the symmetric positive voltage drop as- 
sumed herein, the larger chemical potential is = eV/2. 
This implies that the occupation probability vanishes 
if eV < 2(N — l)U, which is consistent with the occur- 
rence of the steps shown in Fig. [3Jb) . In particular, for 
e\V\ < 2U, only single occupation plays a role. Below we 
will use this fact for establishing an analytical treatment 
in the low-bias regime. 
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Figure 4: (a) Fano factor of the current shown in Fig. 02a) for the 
interaction strength U = 80r. (b,c) Fano factor as a function of the 
detuning £2 = eV ga tc for the bias voltages (b) V = 120r/e and (c) 
V = 280r/e. The dashed line in panel (b) represents the analytical 
result H27H valid in the one-electron regime. 



3.2. Shot noise and Fano factor 

So far we have seen that unless the bias voltage is ex- 
tremely low, the current is always of the order er / fi, i.e. 
changing the gate voltage or the bias voltage modifies the 
current typically by a factor of the order unity. Figure |4] 
demonstrates that the noise, characterized by the Fano 
factor F — iS/e|7|, generally exhibits a more significant 
dependence on both the bias and the gate voltage. 

In the three-electron regime, i.e. for eV > 4U, the Fano 
factor is of order unity and almost independent of e 2 . This 
indicates that the transport process is Poissonian. In the 
two-electron regime, i.e. for U < eV/2 < 2U, the Fano 
factor [Fig. @£c)] exhibits a quadratic dependence on 62 as 
well, but the absolute values are now significantly smaller. 

A rather pronounced dependence on the gate voltage 
Vgatc = £2/e is found in that part of the single-electron 
regime in which the left chemical potential is so large that 
all single particle levels lie within the voltage window, i.e. 
for e 2 < eV/2 < U; see Fig.@Jb). In particular, we observe 
a quadratic growth of the Fano factor, F cx ef, with highly 
super-Poissonian values. This already indicates electron 
bunching, where each bunch consists of roughly F elec- 
trons [40(. 

3.3. Full counting statistics 

A more complete picture can be drawn by considering 
the distribution function of the transported electrons or 
equivalently all cumulants, i.e., the full counting statistics. 
Here we restrict ourselves to the corresponding long-time 
limit which characterizes the low-frequency fluctuations. 
For avalanche- like transport, higher-order cumulants were 
computed 2^, 25, 27, 2^, 23] and recently also measured 
0. 

In a simplified picture, one may consider the electron 
avalanches as particles with charge q = eF that are trans- 
ported in an uncorrelated manner. For this Poisson pro- 
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expressions 



Figure 5: Full counting statistics of the current for t2 = 40r and 
the bias voltages V = 120r/e (one-electron regime) and V = 280r/e 
(two-electron regime) expressed by the cumulants C^. All other pa- 
rameters are as in Fig. [4] Open symbols mark positive values while 
filled symbols correspond to negative values. The line marks the 
behavior for uncorrelated tunneling of charges q = eF in the one- 
electron regime (circles), cf. Eq. H20I I. Inset: Enlargement of the 
lower left corner. 



cess, the cumulants ((5J) grow exponentially fast with their 
index according to the relation [42j 



fe-i 



Ci 



F k ~ l Cx. 



(20) 



Closer theoretical investigations of avalanches with finite 
duration 25j, |3_4|, however, indicate that this picture de- 
serves some refinement. In fact there it was found that 
the cumulants grow even faster. In avalanche diodes, by 
contrast, the opposite was measured, namely that the cu- 
mulants do not grow as fast [4l| . But so far the underlying 
mechanism has not been revealed. 

Figure [5] depicts the super exponential growth of the 
Cfe in both the one-electron regime and the two-electron 
regime. The solid line demonstrates that only the cumu- 
lants of very low order follow the behavior of the Poisson 
process underlying Eq. (I20p . Thus, we also in our case 
observe that already the third cumulant slightly exceeds 
the value P 2 Ci, while higher-order cumulants even assume 
significantly larger absolute values. It is also interesting to 
notice that the sign of the cumulants changes periodically 
with the index k. This hints on the recently conjectured 
"universal cumulant oscillations" 431 ] . 



4. Analytical solution in the one-electron regime 

If the bias voltage is so small that only one electron 
can reside on the triple dot, i.e. for eV < 2U, all relevant 
eigenstates of Ptqd are the empty state |0) and the one- 
electron states. To lowest order in T/e2, the latter are 
given in the basis of the localized states \j) = cj|0) by the 



|0 1 ) = -^(|1)-|3)), 



102) 



1 



\<h) = 



4(T/ £2 )2 
1 



£2 



T 

-1 

£2 



13) 



T 

-3 

£2 



(21) 
(22) 
(23) 



V / l + 2(T/e 2 ) 2 

For symmetry reasons, each of these states couples with 
equal strength to the left and to the right lead, such that 
the corresponding incoherent transition rates from and to 
the leads are given by 



r T 



R,n 



r|«>„|3)| 5 



r, 



(24) 



For the approximate eigenstates (f2"T j) - ([2"3"|) the coupling 
constants read T 1 =T 2 = T/2 and T 3 = rT 2 / e 2- 

If the dot-lead coupling T is sufficiently small, one can 
neglect within a rotating-wave approximation off-diagonal 
elements of the reduced density operator [10( , such that the 
system is well described by the occupation probabilities 



(Po, Pi, P 2 , P 3 ) for the eigenstates |0) and 
The corresponding Liouville operator is 

f-T T x 

ri -r, 
r 2 o 
\r 3 o 



•Cle — 



r 2 
o 

-r 2 
o 



r 3 
o 
o 

-r 3 y 



1,2,3. 



(25) 



where the index "le" refers to the restriction to single oc- 
cupation. We have assumed that the chemical potentials 
are such that all three one-electron states lie within the 
voltage window. Moreover, we have employed the sum 
rule r = J2n=i which follows from the completeness 
relation for the eigenstates \4> n ) in the one-electron sub- 
space. The corresponding particle current operators are 
3~ = 0, while is obtained from C\ e by keeping only 
the elements placed above the diagonal. 

The stationary state of the Liouville operator can now 
be readily computed and reads Po = 1/4 = P„ for all n. 
Inserting this into the current formula (|T71) . yields / = 
eT /Ah. Interestingly enough, the current does not depend 
on the structure of the eigenstates, which is a consequence 
of the mentioned sum rule for the r„. 

For the Liouville operator (j2"5j) . not only the station- 
ary current, but also the zero-frequency noise (I18|) can be 
evaluated exactly. Starting from the latter expression, one 
obtains after some lines of straightforward calculation the 
result 



e 

32h 



3 

E 

n,n' — 1 



r 2 



2 3 

16ft ^ 

n—1 



rir 2 r 3 
r 2 ■ 



(26) 



In the limit Ti. 2 <C T 3 , only terms with T 3 in the denomi- 
nator are relevant, so that we obtain S = (e 2 /ft)r 2 /32r 3 . 
With the above expressions for r„, the corresponding Fano 
factor becomes 



Pie = 



8T 2 



(27) 



It indeed exhibits the predicted parabolic dependence on 
the gate voltage V ga ,tc — ea/e. A quantitative comparison 
of both the current and the Fano factor with the numerical 
results depicted in Fig.Sfb) yields a satisfactory agreement 
in the range considered here. This implies that our one- 
electron model indeed captures the essential features of the 
electron avalanches through the triple quantum dot. More- 
over, it implies that the (average) size of the avalanches is 
q = ee 2 /8T 2 . 

The proportionality of the Fano factor to (e2/T) 2 sug- 
gests for the transport in the one-electron regime the fol- 
lowing scenario: Assume that the off-resonant dot is ini- 
tially unoccupied, while electrons are transported via dots 1 
and 3. Then, according to standard perturbation theory, 
electrons in these dots may tunnel to dot 2 with a probabil- 
ity amplitude proportional to the tunnel matrix element 
divided by the energy difference, i.e. with a probability 
p oc (T/e2) 2 . Thus after on average 1/p electrons have 
been transported, an electron will tunnel to dot 2 and 
cause temporary Coulomb blockade which is resolved only 
when the electron tunnels further. Thus, the dynamics is 
given by periods with an open channel that on average ter- 
minate after an avalanche with 1/p oc (T/e) 2 electrons has 
been transported. This explains the observed proportion- 
ality of the Fano factor F oc (e2/T) 2 . This scenario also 
hints on why the cumulants deviate from conjecture (|20[) : 
The average duration of an avalanche is as long as the 
waiting time between subsequent avalanches, while con- 
jecture (|2T)|) is based on the assumption of a much shorter 
avalanche duration. 

In the limit T3 <C V, it is still possible to evaluate 
some cumulants of higher order analytically, although this 
becomes increasingly tedious. Nevertheless it is worth pro- 
ceeding up to forth order for which we obtain 




Thus, the first two cumulants behave as expected for short 
avalanches with charge q = er/8r3. However, already the 
third and the forth cumulant deviate from the Poissonian 
value by a factor ±3 in compliance with our numerical 
observations in Sec. 13.31 

5. Conclusions 



that is not connected to any lead has the consequence that 
transport through one arm requires co-tunneling. Since 
this reduces the transmission probability of that path, the 
interference pattern as a function of a penetrating flux 
will fade away. We have demonstrated that, nevertheless, 
strongly detuned interferometers bear interesting effects 
manifest in the noise properties of the current. In partic- 
ular, we have shown that strong electron bunching may 
occur. It turned out that this is most pronounced in the 
low bias regime in which Coulomb repulsion forbids the 
occupation of the triple quantum dot by more than one 
electron. 

The prime quantity of interest in that context is the 
Fano factor for which we have predicted huge values: the 
relative noise strength may exceed that of a Poisson pro- 
cess by several orders of magnitude. This has led us to the 
conclusion that the current consists of avalanches with a 
finite duration. The physical reason for this is that elec- 
trons may become trapped in the off-resonant quantum 
dot, such that Coulomb blockade interrupts the trans- 
port until the trapped electron is released. The analysis 
of the higher-order charge fluctuations — the full-counting 
statistics — has revealed that the cumulants grow even su- 
per exponentially with their order. 

For the computation of the full-counting statistics, we 
have employed the iterative scheme recently developed in 
Ref. [33}. This has been essential for the treatment of 
the three-dot problem, since the Hilbert space is already 
too large for a full analytical treatment. Nevertheless, in 
the most relevant regime of single occupation, we have 
performed the first iteration steps analytically, such that 
we have obtained expressions for the first four cumulants 
in the limit of strong detuning. 

In conclusion, we have studied electron interferometers 
in a regime that so far has not attracted much attention, 
most likely due to the lack of pronounced interference ef- 
fects. Precisely in this regime, however, the shot noise 
properties are most interesting and strong electron bunch- 
ing occurs. Thus, quantum dots in ring configuration may 
serve not only for the observation of interference effects, 
but also for the creation of currents with widely tunable 
super-Poissonian fluctuations. 

Acknowledgements 

S.K. likes to thank Peter Hanggi for numerous inspir- 
ing discussions during many years of fruitful collaboration. 
We acknowledge helpful discussions with Christian Flindt. 
This work has been supported by the Spanish Ministry of 
Science and Innovation through project MAT2008-02626, 
a Ramon y Cajal fellowship (S.K.), andaFPI grant (F.D.). 



The emergence of an interference pattern with good 
visibility usually requires the coherent superposition of two 
or more paths that are traversed with like probabilities. 
In an Aharonov-Bohm interferometer formed by quantum 
dots in ring configuration, a significant detuning of the dot 



Appendix A. Rayleigh-Schrodinger perturbation 
theory 

In most textbooks on quantum mechanics, perturba- 
tion theory is formulated for a perturbation linear in the 
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small parameter, while the augmented Liouvillian C x con- 



tains arbitrarily high powers of x; see Eq. (fit))) . In this ap- 
pendix we derive a recursive scheme for the representation 
of an eigenvalue as a Taylor series in a perturbation param- 
eter. Our calculation is inspired by Chap. 5 of Ref. 44j, 
for an alternative derivation see Ref. [45[ . Despite the fact 
that we formulate the problem in terms of a quantum me- 
chanical energy eigenvalue problem, our derivation is not 
restricted to Hermitian operators. 

We consider a "Hamiltonian" H = Ho + V(a) for which 
the perturbation V{a) is an analytical function of the per- 
turbation parameter a and vanishes in the limit a — > 0. 
Thus, it can be decomposed into the series 



V(a) = £Vv fc . 



(A.l) 



k=l 



Our goal is to find a series for the eigenvalue E(a) of H 
that fulfills lim^o E(a) = E , where E is the eigenvalue 
of Hq corresponding to a particular eigenvector \4>o), i.e., 
Hol&o) — E \(j)o)- A central assumption is that both the 
eigenvalue and the corresponding eigenvector |0(a)) can 
be decomposed into a series in a as well, such that we can 
employ the ansatz 



(0o I Q — 0, the relation (0 o |0(a)) = (0o|0o) = 1 holds. An 
important feature of Eq. (|A.5[) is that the second term on 
the right-hand side is of first order in a. Therefore, it can 
be solved iteratively in the following way. 

Multiplying Eq. (|A~4l) by (0 O |, we find that the left- 
hand-side vanishes owing to (0o|-ffo = Eo((f>o |. Thus we 
obtain E(a) — Eq = (4> \V (a)\4>(a)) . We then insert for 
E(a), V(a), and |0(a)) the series (|A.1I) - (|A.3I) and com- 
pare coefficients to obtain the fcth-order energy shift 



(A.6) 



It depends on the still unknown correction of the eigen- 
vector, \4>k-k'), which we determine from Eq. (|A.5|) . We 
again insert Eqs. (IA.1[) - (|A.3[) and compare coefficients to 
find 



\<Pk) 



Q 



H — E[ 



k-l k 

r {j2Ek'\<j>k-k')-Y, v ^-^}- 

k' = l k' = l 



(A.7) 

Equations (jA.6[) and (|A.7|) allow the recursive computation 
of the series (|A.2[) for the eigenvalue shift E(a). 



E(a) =E Q + ) a k E, 



!</>(«)) = I0o) 



fc=i 



a k \(t>k 



(A.2) 



(A.3) 



For the present purpose, it is sufficient to consider only the 
case of non-degenerate Eq. Note that for non-Hermitian 
Hq, (0o I is generally not the Hermitian adjoint of |0o), but 
rather the corresponding left eigenvector, i.e., the solution 
of (0o I #o = -Eo(0o|- We choose the normalization such 
that (0 o |0o> = 1- 

We start by writing the eigenvalue equation H\<j){a)) = 
E(a)\<f)(a)} in the form 

(H - Eo)\4>(a)) = {E(a) - E - F(a)}|0(a)), (A.4) 

which is convenient for the derivation of a formal solu- 
tion and subsequent iteration. Since (Hq — £'o)|0o) = 0, 
Eq. (|A.4|) defines 10(a)) only up to a component propor- 
tional to |0o). Moreover, it implies that the inverse of 
Hq — Eo does not exist. It is however possible to define the 
pseudo-inverse Q(H — E )~ 1 Q [or in short: Q/(H — E )], 
where Q is the projector on the subspace spanned by the 
eigenvectors of Ho with non-zero eigenvalue. In the present 
case of non-degenerate E , it reads Q = 1 — |0 o )(0o|. 

The pseudo-inverse allows one to cast Eq. (|A.4[) into 
the form 



10(a)) = I0o) 



Q 



Ho — Eo 



{E(a)-E -V(a)}\<j>(a)). (A.5) 



The first term on the right-hand side can be multiplied 
by any factor without violating Eq. (|A.4[) . We have cho- 
sen it such that linia^o 10(a)) = |0o)- Moreover, since 
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